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Abstract 

Stable electrostatic levitation and trapping of a neutral, polarizable 
object by a charged ring is shown to be theoretically impossible. Earn- 
shaw's theorem precludes the existence of such a stable, neutral particle 
trap. 

1 Introduction 



In this tribute in honor of the memory of Prof. Dr. Herbert Walther, we consider 
the possibility of extending his famous work on the trapping of an ordered lattice 
of ions [1] in a Paul trap [2] , to the trapping of neutral atoms, and more generally, 
to the possible levitation of a macroscopic neutral polarizable object, in a purely 
electrostatic trap, for example, in the DC electric field configuration of a charged 
ring. Earnshaw's theorem will be extended to the case of such neutral objects, 
and we shall show below that the stable levitation and trapping of a neutral, 
polarizable object, which is a high-field seeker, is generally impossible in an 
arbitrary electrostatic field configuration. We shall do this first for the special 
case of the electrostatic configuration of a simple charged ring, and then for the 
general case of any DC electric field configuration. 

Consider the charged-ring geometry shown in Figure 1. The region near 
the so-called "levitation" point L in this Figure is akin to the focal region of 
a lens in optics. Just as two converging rays of light emerging from a lens 
in physical optics cannot truly cross at a focus, but rather will undergo an 
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Figure 1: A uniformly charged ring with radius a lies on the horizontal x- 
y plane, with its axis of symmetry pointing along the vertical z axis. Can 
levitation and trapping of a neutral particle occur stably near point L, where 
there is a convergence of E-ficld lines? 
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"avoided crossing" near the focal point of this lens due to diffraction, so likewise 
two converging lines of the electric field cannot cross, and therefore they will 
also undergo an "avoided crossing" near L. There results a maximum in the 
z component of the electric field along the vertical z axis at point L. The 
resulting "avoided crossing" region of electric field lines in the vicinity of point 
L is therefore similar to the Gaussian beam-waist region of a focused laser beam. 
Ashkin and his colleagues [3] showed that small dielectric particles, which are 
high-field seekers, are attracted to, and can be stably trapped at, such Gaussian 
beam waists in "optical tweezers" . Similarly here a neutral dielectric particle, 
which is a high-field seeker, will also be attracted to the region of the convergence 
of E-ficld lines in the neighborhood of L, where there is a local maximum in the 
electric field along the z axis. The question arises: Can such a high-field seeker 
be stably levitated and trapped near L? 



2 Calculation of the electric potential and field 
of a charged ring 

The electric potential at the field point P due to a charge element dq' of the 
ring is given in by 

d$ = # (1) 
r 

where the distance r from the source point, whose coordinates are (x', y' , 0), to 
the field point P, whose coordinates are (x,0, z), is 



r = ^J{x' - x) 2 + y' 2 + z 2 . (2) 

(Primed quantities refer to the source point; unprimed ones to the field point). 
Since the charged ring forms a circle of radius a which lies on the horizontal x-y 
plane, 

x ' 2 + y ' 2 = a 2 . (3) 

An infinitesimal charge element dq 1 spanning an infinitesimal azimuthal angle 
of d<j) can be expressed as follows: 

dq' = (*-) a # ' = - (*-) ad{x ' /a) (4) 
\2iraJ v \2TTaJ ^1 - ( x '/ a ) 2 V ' 

where Q is the total charge of the ring. Let us introduce the dimensionlcss 
variables / 

S = — > V = — . C = e = -• 5) 
a a a a 

Thus 

Q 

dq 



2 71 



3 



Due to the bilateral symmetry of the ring under the reflection y' — ► —y', it is 
useful to sum up in pairs the contribution to the electric potential from symmet- 
ric pairs of charge elements, such as dq[ and dq' 2 with coordinates (x',+y',0) 
and (x',—y',0), respectively, shown in Figure 1. These two charge elements 
contribute equally to the electric potential $ if they span the same infinitesimal 
azimuthal angle d(f>' . Thus one obtains 

$ (e, C) = — / +1 dg ; 1 — r 1 . (6) 

Along the z axis, this reduces to the well-known result 

$( £ = 0,0 = --7=Uj = -p£=j ■ (7) 

a VT + T V^ 5 + a 2 

The z component of the electric field, which is the dominant .E-field component 
in the neighborhood of point L, is given by 

9$ Q 1 1 



9z 



C / # ^ =, — . (8) 



-i Vi-e' 2 ( v /e2_ 2e e' + l + C 

Along the z axis, this also reduces to the well-known result 

E z = QZ , (9) 

(z 2 + a 2 f /2 

which has a maximum value at 

Zo = 7! orCo = 7r (10) 

The "levitation" point L then has the coordinates 

l(o,o,^), (11) 

neglecting for the moment the downwards displacement of a light particle due 
to gravity. 

The potential energy U for trapping a neutral particle with polarizability a 
in the presence of an electric field (E x , E yi E z ) is given by 



U = -\a(El+El+Et) ^-\aE 2 z , 



since the contributions to U from the x and y components of the electric field, 
which vanish as e near the z axis for small e, can be neglected in a small 
neighborhood of L. 



4 



We now calculate the curvature at the bottom of the potential-energy well 
U along the longitudinal z axis, and also along the transverse x axis. The force 
on the particle is given by 

F = -VU . (13) 
Therefore the z component of the force is, to a good approximation, 

dE 

F z = aE z —^, (14) 
oz 

and the Hooke's law constant k z in the longitudinal z direction is given by 

0F Z ffdE z \ 2 d 2 E z \ 

where all quantities are to be evaluated at L where e = and £ = \j\f2. 
Taking the indicated derivatives and evaluating them at L, one obtains 

7, I 32 aQ 2 

k ^ = + 81— ' (16) 
where the positive sign indicates a longitudinal stability of the trap in the ver- 
tical z direction. 

The x component of the force is, to the same approximation, 

F. = aE,°£, (17) 
and the Hooke's law constant k x in the transverse x direction is 

3F X \ fdE z \ 2 „ d 2 E z \ 

*- = -ftf = -°{Ur) +E '-toF] ' (18) 

where again all quantities are to be evaluated at L where e = and Co = V \/2- 
Again taking the indicated derivatives and evaluating them at L, one obtains 

M, = --^ (19) 

x\l 81 a 6 ' y ' 

where the negative sign indicates a transverse instability in the horizontal x 
direction. 

Similarly, the Hooke's law constant k y in the transverse y direction is 

16 aQ 2 
81 

where the negative sign indicates a transverse instability in the horizontal y 
direction. Note that the trap is azimuthally symmetric around the vertical axis, 
so that the x and y directions are equivalent to each other. Because of the 
negativity of two of the three Hooke's constants k x , ky, and k z , the trap will 
be unstable for small displacements in two of the three spatial dimensions near 
L, and hence L is a saddle point. Note also that the sum of the three Hooke's 
constants in Equations (16), (19), and (20) is zero, i.e., 

k X + ky+ k, = 0. (21) 
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3 Earnshaw's theorem revisited 



We shall see that Equation (21) can be derived from Earnshaw's theorem when 
one generalizes this theorem from the case of a charged particle to the case of a 
neutral, polarizable particle in an arbitrary DC electrostatic field configuration. 
A quantitative consideration of the force on the particle due to the uniform 
gravitational field of the Earth, in conjunction with the force due to the DC 
electrostatic field configuration, does not change the general conclusion that 
the mechanical equilibrium for both charged and neutral polarizable particles is 
unstable. 



4 Charged particle case 



We shall first briefly review here Earnshaw's theorem [4], which implies an 
instability of a charged particle placed into any configuration of electrostatic 
fields in a charge-free region of space in the absence of gravity. Suppose that 
there exist a point L of mechanical equilibrium of a charged particle with charge 
q in the presence of arbitrary DC electrostatic fields in empty space. The 
potential $ for these fields obey Laplace's equation 



<9 2 $ 



<9 2 $ 
dy 2 



<9 2 $ 



dx 2 dy 2 dz 2 
Now the force on the charged particle is given by 



= 0. 



<9$ 

' dx 
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(F x ,F y ,F z ) 



(22) 



(23) 



where e x , e y , e z are the three unit vectors in the x, y , and z directions, respec- 
tively. By hypothesis, at the point L of mechanical equilibrium 







<9$ 


d<5> 


dx 
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dy 


L dz 



= 



(24) 



Stable equilibrium would require all three Hooke's constants k x , k y , and k z at 
point L to be positive definite, i.e., 



2? n 



dF x 
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dx 2 
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dF z 
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L 



> 



> 



> 0. 



(25) 
(26) 
(27) 
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However, Laplace's equation, Equation (22), can be rewritten as follows: 

k x + k y + k z = 0, (28) 

i.e., the sum of the three components of Hooke's constants for the charged par- 
ticle must be exactly zero. The simultaneous positivity of all three Hooke's 
constants is inconsistent with this, and hence at least one of the Hooke's con- 
stants along one of the three spatial directions must be negative. Therefore the 
system is unstable. 



The azimuthally symmetric field configurations like that of a charged ring 
is an important special case. Let z be the vertical symmetry axis of the ring. 
Suppose that there is stability in the longitudinal z direction (such as along the 
z axis above point L), so that 

k z > . (29) 

By symmetry 

k x = k y = k± (30) 
so that Equation (28) implies that 

k± = -h z < , (31) 

implying instability in the two transverse x and y directions. 

Conversely, suppose there is instability in the longitudinal z direction (such 
as along the z axis below point L), so that 

k z < . (32) 

Again, by symmetry 

= ky — kj_ (33) 

so that Equation (28) implies that 

k± = ~k z > , (34) 
implying stability in the two transverse x and y directions. 



5 Adding a uniform gravitational field such as 
the Earth's, in the case of a charged object 

The potential energy of a charged, massive particle in a DC electrostatic field 
in the presence of Earth's gravitational field is 

Utot = q$ + mgz. (35) 
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Note that the term due to gravity, i.e., the mgz term, is linear in z, and therefore 
will vanish upon taking the second partial derivatives of this term. Therefore 
the Hooke's constants k x , k y , and k z will be unaffected by Earth's gravity. The 
force on the particle is 

Ftot = -Wtot = -<?V$ - mge z (36) 

where e z is the unit vector in the vertical z direction. In equilibrium, F tot = 
0, but this equilibrium is again unstable, since upon taking another partial 
derivative of the term mge z with respect to z will yield zero, and therefore 
all of the above Hooke's law constants are the same in the presence as in the 
absence of Earth's gravity. 



6 Generalization to the case of a neutral, polar- 
izable particle 



Now suppose that there exists a point L of mechanical equilibrium of the neu- 
tral particle with positive polarizability a > somewhere within an arbitrary 
electrostatic field configuration. Such a particle is a high-field seeker, and hence 
point L must be a point of high field strength. Choose the coordinate system 
so that the z axis is aligned with respect to the local dominant electric field 
at point L. Thus the dominant electric field component at L is thus E z . The 
potential energy U for a neutral particle with polarizability a in the presence 
of an electric field (E x , E y ,E z ) is given by 



U 



1 , 9 



Ez 



1 9 

2 2 



(37) 



since the contributions to U from the x and y components of the electric field, 
which vanish as e near the z axis for small s, can be neglected in a small 
neighborhood of L. The force on the particle is 



F = -VU . 

Therefore the z component of the force is, to a good approximation, 

dE z 



F z = aE z 



dz 



and the Hooke's law constant k z in the z direction is given by 



dF z 
dz 



(dE z 



\ dz 



E z 



b 2 e z 

dz 2 



= -a E. 



d 2 E z 



dz 2 



(38) 



(39) 



(40) 



where the last equality follows from the hypothesis of mechanical equilibrium 
at point L. 



Similarly, the x component of the force is, to the same approximation, 

dE z 



F x = aE z 



dx 



and the Hooke's law constant k x in the x direction is given by 



k x = = -a 

ox 



dEz 
dx 



E z 



d 2 E z 
dx 2 



= -a E z 



o 2 e z 



dx 2 



(41) 



(42) 



where the last equality follows from the hypothesis of mechanical equilibrium 
at point L. 

Similarly the y component of the force is, to a good approximation, 

dE z 



F y = aE z r 
dy 

and the Hooke's law constant k y in the y direction is given by 



dF z 
dy 



(dE z 



E, 



d 2 E z 
dy 2 



-a E. 



d 2 E z 



dy 2 



(43) 



(44) 



where again the last equality follows from the hypothesis of mechanical equilib- 
rium at point L. 

Thus the sum of the Hooke's law constants along the x, y, and z axes is 
given by 

'd 2 E z d 2 E z d 2 E z 



-a j£ 

-a E 



dx 2 dy 2 dz 2 
dz 2 



Therefore 



dz { dx 2 dy 2 



k„ + K)\ T =0 , 



(45) 



(46) 



(k x I ny I n Z )\ L 

and again, the sum of the three Hooke's law constants must be exactly zero 
according to Laplace's equation. 

Suppose that the system possesses axial symmetry around the z axis with 



k z > , 

i.e., with stability along the z axis. Then by symmetry 



ky — k_[ 



so that Equation (46) implies that 



k± = --k z < 



(47) 
(48) 

(49) 



implying instability in both x and y directions. This is exactly what we found 
by explicit calculation for the case of a neutral, polarizable object near point L 
of the charged ring. 
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7 Adding a uniform gravitational field such as 
the Earth's, in the case of a neutral, polariz- 
able object 



The potential energy of a neutral, polarizable, massive particle in a DC electro- 
static field plus Earth's gravity is 

Utot = U + mgz. (50) 

Again, note that the term due to gravity, i.e., the mgz term, is linear in z, and 
therefore will vanish upon taking the second partial derivatives of this term. 
Therefore again the Hooke's constants k x , k y , and k z will not be affected by 
Earth's gravity. The force on the particle is 

Ftot = -VUtot = -qVU - mge z (51) 

where e 2 is the unit vector in the vertical z direction. In equilibrium, F to t = 
0, but this equilibrium is again unstable, since upon taking another partial 
derivative of the term mge z with respect to z will yield zero, and therefore 
again all of the above Hooke's law constants are the same in the presence as in 
the absence of Earth's gravity. 

8 Ways to evade Earnshaw's theorem 

Some known ways to evade Earnshaw's theorem and thereby to construct a truly 
stable trap for charged or for neutral particles are (1) to use non-electrostatic 
fields such as a DC magnetic field (e.g., the Penning trap [5]) in conjunction 
with DC electric fields, or (2) to use time-varying, AC electric fields, rather 
than DC fields (e.g., the Paul trap [2]), or (3) to use active feedback to stabilize 
the neutral equilibrium of a charged particle in a uniform electric field, such 
as was done for a charged superfluid helium drop [6], or (4) to use the low- 
field seeking property of neutral, diamagnetic objects to levitate them in strong, 
inhomogencous magnetic fields [7] . The latter two methods may be practical for 
levitating the superfluid helium "Millikan oil drops" in the experiment described 
in [8]. 
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